Scattering phase of quantum dots: Emergence of universal behavior 
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We investigate scattering through chaotic ballistic quantum dots in the Coulomb blockade regime. Focusing 
on the scattering phase, we show that large universal sequences emerge in the short wavelength limit, where 
phase lapses of n systematically occur between two consecutive resonances. Our results are corroborated by 
numerics and are in qualitative agreement with existing experiments. 
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Quantum mechanics fundamentally differs from classical 
mechanics in that time evolutions are determined by com- 
plex probability amplitudes instead of real probabilities. The 
associated phase is a key element to understand mesoscopic 
transport experiments on Aharonov-Bohm (AB) conductance 
oscillations, weak localization, and conductance fluctuations 
However, these transport measurements, like any other 
measurement, do not directly measure scattering phases. In 
their phase-sensitive experiment, Yacoby et al. pioneered an 
entirely new field of mesoscopic physics, by embedding a bal- 
listic quantum dot (QD) in one arm of an AB interferometer 
J2l- The sustained interest in these and following lUHt] ex- 
periments, that persists until today, relies on the difficulty to 
consistently and generically explain these measurements. 

When a QD in the Coulomb-blockade (CB) regime is 
placed in one arm of a mesoscopic ring threaded by a flux 
4> (see the inset of Fig. Q} the conductance through the ring 
reads mlE 



(1) 



with the quantum of flux 0o- Measuring the lowest harmon- 
ics of the AB oscillations in <p allows to extract the phase f}\ 
as a function of a gate voltage V g applied to the QD. Under 
suitable conditions f}\ can be related to the scattering phase 
of the QD 17]. A multi-terminal configuration is required in 
order to get a one-to-one relationship between the two phases, 
while in a two-terminal set-up j3\ can only take the values 
or n In the multi-terminal case a gradual increase of n 

in the transmission phase is obtained as a function of the gate 
voltage for every CB peak, in agreement with the Friedel sum 
rule H8l— [ 1 Oil . Abrupt lapses of n occur between resonances in 
both configurations at values of V g for which the transmis- 
sion amplitude of the QD is so small that the AB oscillations 
are below the experimental visibility threshold. For relatively 
large QD with a few hundred electrons |2j,|3[], the lapses were 
seen to appear systematically between every consecutive pair 
of resonances. This behavior, termed universal, is observed 
in two- or multi-terminal configurations. To the best of our 



knowledge, it is still unexplained. Smaller dots were more re- 
cently investigated, where the number of electrons was tuned 
from about twenty down to zero |5[ . A crossover from the uni- 
versal regime, where successive peaks are in phase, to a meso- 
scopic regime, where phase lapses occur in a random fashion, 
was observed when decreasing the number of electrons in the 
QD by the action of Vg. Multi-terminal configurations with 
one QD in each arm of the AB interferometer have also been 
investigated [4], which determined the important role of the 
magnetic field in the phase lapses. The understanding of the 
crossover from the mesoscopic to universal regime and the 
role of the symmetries are the main goals of our work. 

For lateral QD in the CB regime a single lead channel dom- 
inantly couples to the QD, and transport through the QD is 
characterized by a 2 x 2 scattering matrix 



S = 



r f 
t r' 



ie^ cos 9 e' v sin 6 
e~ lT> sin 6 ie~% cos 6 



(2) 



We note t (f) and r (/) the transmission and reflection ampli- 
tudes, respectively, for particles coming from the left (right) 
of the QD. The angle 9 is restricted to the interval [0, n), while 
the phases a, r\ and £ are defined on [0, 27r). The scattering 
phase a is related to the density of states of the QD through 
the Friedel sum rule |0, fsl— T]~qT| . When considering the phase 
evolution as a function of an external parameter (like V g ) it 
is convenient to work with the accumulated phase a c , whose 
range of definition is not restricted to the interval [0, 2n). 

Eq. (f2]) represents the most general 2x2 unitary matrix. 
Time-reversal invariance implies that 77 = or n, while the 
right-left parity symmetry fixes £ to either or n. When 
the many-body problem is considered in its full complexity, 
S represents an effective one-particle scattering matrix that 
can be obtained for instance through the embedding method 
I U , 12] • If we restrict the description of the many-body prob- 
lem to the constant-interaction model (CIM) the CB phenom- 
ena can be interpreted in terms of single-particle quantities 
1 1311 . and close to resonances, S is given by the transmission 
and reflection amplitudes of single-particle states in a mean- 
field potential. When the underlying classical scattering is 



2 



chaotic, S has well defined statistical properties determined 
by the symmetries of the problem only [14]. In structures with 
time-reversal symmetry, t = e' a sin 9 and, even if t is a com- 
plex continuous function of the real variable V g , the phase a 
exhibits a phase lapse of n whenever t vanishes. 

Within the CIM and in the absence of a magnetic field the 
one-particle wave-functions can be chosen to be real. With the 
exception of peculiar cases with very different values of the 
resonance widths between consecutive resonances, a zero of 
the transmission generically appears between two successive 
resonances depending on the sign of II 1 C)fl 



D n = TiYnTn+lfn+l 



(3) 



The partial-width amplitude (or effective coupling) of the 
eigenstate i/r„ of the QD is given by 11311 



i n 



H 2 kP c 



1/2 



Jo 



dy® m (y)Mx m ,y)- (4) 



Here, <I>i( r ) is the first transversal subband wave-function in 
the the left (right) lead of width W, the integration is along 
the transverse coordinate y at the entrance or exit of the QD 
located at x — x 1 ®, P c is the transparency of the tunnel bar- 
riers, k is the Fermi wave-vector in the leads, and m is the 
electron mass. In lattice models with one-dimensional leads 
10, OH the two partial-width amplitudes are simply propor- 
tional to the value of the wave-function at the extreme point 
connecting the QD to the corresponding lead. 

We call y\y x n the parity of the n th resonance lflQIl . When 
D n > (equal parity of the n th and (« + l) st resonances) 
there is one zero (or an odd number of zeros) between the 
n th and (n + l) st resonances, while for D„ < (opposite par- 
ities) there is no zero (or an even number of zeros) between 
the two resonances. For one-dimensional scatterers, D n is al- 
ways negative, which results in the absence of transmission 
zeros reflecting the impossibility of obtaining destructive in- 
terfering paths in one dimension|8|]. Numerical simulations 
on two-dimensional disordered lattice systems, on the other 
hand yield an equal probability for positive or negative D„ 
ifioll . which is at odds with the experimental observation of a 
universal regime with zeros between any two consecutive res- 
onances. Several theoretical refinements have been proposed 
to solve the puzzle lflol[l6l -23l. where specific geometries or 
effective couplings, or extensions of the CIM are considered. 

The latter path is in principle the most natural one, and was 
actually already suggested in Ref. [2]. However, the small 
sizes that can be handled within a full many-body descrip- 
tion make the universal regime hardly reachable. Moreover, 
even if the influence of electron-electron interaction can be 
spotted in some circumstances the correspond- 

ing results are not generic. On the other hand, the CIM gives 
an excellent description of the statistical distribution of the 
height of the CB peaks (l3, 15], which is determined by the 
one-particle widths y]f. It is then expected that the statistics 
of the resonance parities and consequently of the transmis- 
sion zeros are within the reach of the CIM. Below, assuming 



0.5 




-0.5 



FIG. 1: th order Bessel function J (kL). The shaded areas indi- 
cate regions where D„ is negative. For two-dimensional systems, the 
width of these intervals, and the probability of not obtaining a phase 
lapse between two resonances, decrease as LAk ~ (kL)~ l . Inset: 
Aharonov-Bohm interferometer, threaded by a flux <p, with an asym- 
metric dot (shaded) tunnel-embedded in its upper arm. The distance 
between the entrance and exit points of the dot is L. The dashed lines 
on the arms of the interferometer stand for any number of potential 
additional leads. 



that this is the case, we incorporate well known correlations 
of quantum chaotic wavefunctions in our analysis of D„. In 
this way we account for the transition from a mesoscopic to a 
universal regime. 

Long-range wavefunction correlations in quantum chaotic 
systems were first pointed out by Berry [24], who suggested 
to model the wavefunctions as random superpositions of plane 
waves with fixed energy ft 2 k 2 /(2m). For a two-dimensional 
chaotic billiard with eigenfunctions i//„, this gives the wave- 
function correlator 



^„(r)^„(r') = J (k\r - r'\)/& . 



(5) 



The bar stands for a local average, tft is the area of the billiard 
and Jo is the th Bessel function of the first kind. Corrections 
to Eq. (0 appear in confined systems when the observation 
points approach the boundaries [25], which is the case we are 
interested in. For r + r', these corrections are however small 
and we will neglect them. Eq. (0 has been successfully used 
to explain the long-range (in energy) modulation of the peak- 
height distribution in the CB regime 12^1 . 

If we call L the distance from the entrance to the exit of 
the QD and assume that successive eigenfunctions are uncor- 
rected we have for the case W <k L of relevance for us 



D„ ~ J (k„L)J (k n+l L) . 



(6) 



Taking L as the typical linear dimension of the QD, we have 
K+\ - k„ ^ Ak = nj(kl}). Thus, D„ is negative when the 
Bessel functions in Eq. © have different signs, which only 
happens if kL falls in an interval of length LAk before one 
zero of Jo(kL). This is sketched in Fig.Q] In the semiclassical 
limit kL » 1 we have Jo(kL) ~ V2 / (nkL) cos (kL - n/4) and 
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FIG. 2: Transmission \t\ 2 and scattering phase a as a function of the 
applied gate voltage V g . The number of electrons in the dot varies 
between 174 and 182. Smooth increases of n in a obtained for each 
resonance alternate with phase lapses of —n when the transmission 
vanishes (dotted vertical lines). 



therefore the probability f of obtaining a negative D n can be 
estimated as the ratio of 2LAk over the period In, that is, 



kL 



(7) 



This simple analysis explains why when the dot is progres- 
sively filled the equal-parity case (D„ > 0) results with prob- 
ability approaching one. Moreover, it also shows that the se- 
quence of peaks and zeros appears over intervals in kL that 
are of length n, containing a number kL of resonances. The 
sole assumption that wave-functions have quantum chaotic 
correlations thus predicts the emergence of large universal se- 
quences of resonances and transmission zeros. This is our 
main result. 

In order to illustrate and test our analytical prediction we 
have performed numerical calculations of the transmission 
amplitude of spinless electrons through a noninteracting QD 
with the shape of a desymmetrized stadium billiard, where 
one of the quarter circles is replaced by a cosine curve. This 
is sketched in the inset of Fig. Q] The QD is connected to 
leads through tunnel barriers as in Ref. 11511 . The parame- 
ter of variation is the gate voltage V g which induces changes 
in the wave-vector k within the QD. Phase increases of n are 
correlated with transmission peaks, while phase lapses of n 
are assigned at the values of V g where the transmission van- 
ishes. We show in Fig. [2] a sequence of eight resonances for 
174 to 182 electrons on the dot, similar to the experiments of 
Refs. HH. This sequence exhibits a perfect alternation of 
resonances and zeros. We found that this is the generic behav- 
ior for kL > 10. 

We next present in Fig. [3] (thick lines) the accumulated 
scattering phase (the superscript indicates that we fol- 
low the convention used in the experiments of taking all phase 
lapses as —jf), together with the number of accumulated res- 
onances (zeros) N r ( Z ). N t and N z grow with almost the same 
mean rate for kL > 10 (at least up to the maximum values of 
kL =! 100 that we used). In the inset we show the probabil- 
ity V = (AN r - AN z )/AN r of obtaining more resonances than 
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FIG. 3: Thick lines: accumulated scattering phase a^/n, and num- 
ber of resonances (zeros) iV r(z ) as a function of V g (or kL). The shaded 
spots indicate the regions where N z lags N t . They are separated ap- 
proximately by 7r in kL. Thin lines: a c /n for a small positive and 
negative magnetic field. Inset: P = (AiV r - AAQ/AiVr taken at vari- 
ous intervals of V g as a function of kL. The line is a guide-to-the-eye 
that decreases as (kLy 1 showing the good agreement of the numeri- 
cal calculations with the prediction of Eq. 10. 



transmission zeros in an interval of Vg. f is proportional to 
the difference of slopes AN T /AV g - AN z /AV g . In agreement 
with the theoretical analysis, the regions where there are no 
zeros between resonances show up with a periodicity (in kL) 
of approximately ji and are separated by intervals where there 
is an alternation of peaks and zeros. Upon increasing k, se- 
quences exhibiting perfectly alternating resonances and zeros 
become larger and larger, favoring the observation of the uni- 
versal behavior. 

The importance of a magnetic field in the phase-sensitive 
measurements has been underlined in experimental (4)] and 
theoretical |8, 27] works. Once we have a magnetic field 
breaking time-reversal symmetry, the wave-functions are no 
longer real, D„ becomes complex, and our analysis in terms 
of the parity of the resonances is no longer applicable. We in- 
cluded a magnetic field in our numerical calculations, using a 
desymmetrized structure in order to truly break time-reversal 
symmetry M28I1 . We observed that the exact transmission zeros 
obtained at B = become minima with small (but finite) val- 
ues. Such a behavior is expected, since a mono-channel QD 
has divergent probability of exhibiting a vanishing \t\ when 
B = 0, while the transmission distribution is uniform in the 
interval (0, 1) for the case where time-reversal symmetry is 
completely broken II 1 31 - 4 1 511 . 

Avoiding transmission zeros in the complex plane elimi- 
nates the n phase lapses in favor of continuous jumps with 
large (but finite) derivatives and magnitude < n. At finite 
fields the ambiguity in the definition of the accumulated phase 
a c is then lifted. In Fig. [3] we present a c for a very small pos- 



4 



itive (negative) field (thin lines). The phase jumps for small 
positive fields are opposite to those for negative fields and the 
difference in the accumulated phases is of statistical nature (of 
the order of the square root of the number of avoided zeros), 
illustrating the diffusion away from the origin in the complex 
plane for small B fields. On the other hand, the mean slope 
of the accumulated phase is the same in both cases, and coin- 
cides with the slope of jrN r . We point out that the difference 
between the Friedel and transmission phases 1 8, 9] arises from 
the ambiguity in the definition of a c at B — 0. This ambiguity 
is lifted if we define a c (B — 0) = \m\B^,o*a c {B). 

Experimentally, the visibility threshold for the AB oscilla- 
tions can make the continuous phase evolution at small fields 
indistinguishable from the phase lapses at B = 0. Breaking 
time-reversal symmetry reduces the probability of obtaining 
very small transmission values and thereby favors the observa- 
tion of a continuous phase evolution as a function of V g . The 
magnetic field needed for breaking time-reversal symmetry in 
the QD scales as k' 1 ? 2 JT3L fl3f1 . Experimentally, larger fields 
have been used in Ref. [2], where however abrupt phase lapses 
n are enforced by the use of a two-terminal setup. Refs. i3l|5[] 
on the other hand used fields too weak to break time-reversal 
symmetry. Ref. [4] reported phase lapses, and thus univer- 
sal behavior, in specific magnetic field ranges only, where, 
presumably, the transmission drops below the experimental 
visibility threshold and is thus indistinguishable from a true 
zero. We predict that, under the experimental conditions of 
Refs. [3J,|5[] a magnetic field of the order of 500 G will elimi- 
nate some of the phase lapses observed at weak fields. 

Another test would be to investigate diffusive QD. In a 
two-dimensional system with short-range disorder character- 
ized by an elastic mean-free-path t <K L, the right-hand side 
of Eq. (0 is suppressed by an exponentially small prefactor 
e~ L / 2e , When there is a sufficient amount of disorder, one ob- 
tains V — 1/2, as in Ref. [10], and thus the universal regime 
disappears. The experiments report values I ^ 5 — 15/rni yl- 
5] with QD's of submicron sizes, well in the ballistic regime. 
We predict that a larger, more disordered QD would exhibit 
a mesoscopic regime with phase lapses randomly distributed 
between consecutive resonances with a probability f = 1/2. 
We also note that the weaker screening of impurities in few- 
electron QD makes the dots less ballistic, favoring the obser- 
vation of the mesoscopic regime. 

Up to now we have ignored the effects of the electronic spin. 
The first such effect is a trivial factor of 2 in the density of 
states that affects the relation between k and N. Second, within 
the CIM the filling of spin-degenerate one-particle states is by 
pairs lfl3ll . We note that this spin degeneracy would divide by 
two the probability of obtaining, in our analytical and numeri- 
cal approaches, consecutive resonances without a phase lapse, 
and hence favor the universal behavior. 

Our single-particle theory for ballistic QD thus predicts the 
emergence of a universal behavior of scattering phases at large 
kL. Our numerics indicate that this crossover occurs around 
kL =t 10 which, for two-dimensional structures, correspond 
to putting ~ 15 electrons on the dot. This is qualitatively in 



agreement with the crossover reported in Ref. and with the 
universal behavior reported in Refs. 

fid 

, which work in the 

range kL > 50. We note that another condition is that kL is 
large enough that the wavefunctions resolve the chaoticity of 
the cavity, which also usually occurs around kL 10. We 
predict the disappearance of the observed universal regime (i) 
in the presence of a large magnetic field and (ii) in larger, more 
disordered dots in the diffusive regime. 

In conclusion we have provided quantitative, checkable 
predictions for the probability of observing long sequences 
of alternating transmission zeros and resonances in scattering 
through quantum dots, which are consistent with the experi- 
ments of Refs. IH-Hl- We stress the probabilistic character of 
our findings, and in particular that the absence of phase lapses 
between resonances is always possible. We hope that this will 
stimulate new experimental investigations. 
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